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Abstract. The Calogero-Sutherland model occurs in a large number of physical contexts, 
either directly or via its eigenfunctions, the Jack polynomials. The supersymmetric coun¬ 
terpart of this model, although much less ubiquitous, has an equally rich structure. In 
particular, its eigenfunctions, the Jack superpolynomials, appear to share the very same 
remarkable combinatorial and structural properties as their non-supersymmetric version. 
These super-functions are parametrized by superpartitions with fixed bosonic and fermionic 
degrees. Now, a truly amazing feature pops out when the fermionic degree is sufficiently 
large: the Jack superpolynomials stabilize and factorize. Their stability is with respect to 
their expansion in terms of an elementary basis where, in the stable sector, the expansion 
coefficients become independent of the fermionic degree. Their factorization is seen when 
the fermionic variables are stripped off in a suitable way which results in a product of two 
ordinary Jack polynomials (somewhat modified by plethystic transformations), dubbed the 
double Jack polynomials. Here, in addition to spelling out these results, which were first 
obtained in the context of Macdonal superpolynomials, we provide a heuristic derivation of 
the Jack superpolynomial case by performing simple manipulations on the supersymmetric 
eigen-operators, rendering them independent of the number of particles and of the fermionic 
degree. In addition, we work out the expression of the Hamiltonian which characterizes 
the double Jacks. This Hamiltonian, which defines a new integrable system, involves not 
only the expected Calogero-Sutherland pieces but also combinations of the generators of an 
underlying affine SI 2 algebra. 

Key words: Jack polynomials; supersymmetry; Calogero-Sutherland model; integrable 
quantum many-body problem; affine algebra 

2010 Mathematics Subject Classification: 05E05; 81Q60; 81R12; 37J35 


1 Introduction 

The Calogero-Sutherland (CS) model is an integrable quantum many-body problem with a sur¬ 
prisingly vast number of connections with different physical problems^. Restricting ourself to 
the last 15 years, we can identify a number of important applications of the model itself or its 
eigenfunctions, the Jack polynomials. These are symmetric polynomials labeled by a partition 
A = (Ai,..., Aat) and depending upon a free parameter a (which is the inverse of the CS coupling 
constant). 

*This paper is a contribution to the Special Issue on Exact Solvability and Symmetry Avatars in honour of 
Luc Vinet. The full collection is available at http://www.emis.de/journals/SIGMA/ESSA2014.html 

^For a representative sample of such connections and applications found in the period 1970-1995, see [45]. 
A short historical sketch of the fluctuations in the interest for these models, up to circa 2000, is presented in the 
introduction of [13]. 
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Early in the targeted period, a remarkable property of these polynomials has been unraveled 
by mathematicians [20]. For special partitions, namely with parts satisfying Aj —Aj+fc > r and for 
a = —{k + \)/{r — 1), the Jack polynomials vanish whenever A; + 1 of their variables are equal. 
This technical (dubbed clustering) property made them useful as trial wavefunctions for the 
fractional Hall effect [8]. In this context, the restriction on sequences of A: + 1 contiguous quasi¬ 
particle modes (the parts of the partition) is viewed as a generalized version of the Pauli exclusion 
principle. In conformal field theory (CFT), these special Jacks correspond to the polynomial 
part of the correlators of N generalized parafermionic fields associated to the IF(A; + 1, A; -|- r) 
minimal models [7, 19]. 

A completely different connection with correlators in CFT is displayed in [11, 17]. Still in 
the context of CFT, the remarkable connection between Virasoro singular vectors and Jack 
polynomials found in [4, 5, 28] has been further clarified recently in [34, 36]. The technology 
of Jack polynomials can even be used to derive the spectrum of the Virasoro minimal mo¬ 
dels [34, 46]. These applications have recently been lifted to the s[(2) WZW model at fractional 
level [35]. 

But perhaps the most spectacular application of the Jack polynomials in CFT is linked 
to the Alday-Gaiotto-Tachikawa (ACT) correspondence [3]. This relates conformal blocks in 
two-dimensions CFT to the U(2) Nekrasov instanton partition function [31] pertaining to four¬ 
dimensional SU(2) supersymmetric Yang-Mills theory. The Jack polynomials have been shown 
to be key components of a new AGT-motivated basis of states in 2d-GFT [2]. More precisely, the 
Jack polynomials appear there in a generalized version which is indexed by a pair of partitions 
and decomposes into product of two Jacks with different arguments [2, 30]. 

Here we present a somewhat analogous type of generalization of the Jack polynomials also 
labelled by two partitions. These new generalized Jacks arise directly from the construction 
of the supersymmetric counterparts of the Jack polynomials, the Jack superpolynomials [14]^. 
The latter are eigenfunctions of the supersymmetric version of the GS model [40]^. It tuns out 
that for excited states with large fermionic degree, the eigenfunctions acquire an unexpected 
stability behavior. More remarkably, in this stability sector, these eigenfunctions (after a minor 
transformation) factorize into a product of two Jack polynomials. This factorization is highly 
non-trivial: there is a sort of twisting in the coupling constant (the free parameter a), which 
is different for the two constituent Jacks, and a reorganization of the variables (technically: 
a plethystic transformation)'^. The factorized form of the eigenfunctions is referred to as the 
“double Jack polynomials”. We stress that the non-trivial structure of these double Jacks is 
inherited from the supersymmetric construction, which thus serves as a bridge linking the Jacks 
to their double version. 

These peculiar properties of stability and factorization have first been observed at the level 
of the Macdonald generalization of the Jack superpolynomials [9]. Here we make explicit the 
one-parameter limit characterizing the Jacks. However, instead of adapting to the Jack case the 
rather formal and technical argument of [9], we provide a simple heuristic argument leading to 
the stability and factorization properties. Basically, the strategy amounts to performing simple 
transformations of the two defining eigen-operators of the supersymmetric GS model in order 
to make them independent of the number of variables and the fermionic degree. In addition, we 
unravel the integrable structure underlying the double-Jacks by constructing the Hamiltonian 


^We also use the terminology “Jack polynomials in superspace”. For the motivation underlying this choice 
and the (non-)relationship of these objects with similarly-named polynomials discussed in the literature, see 
Appendix A. 

^The Jack superpolynomials also have clustering properties similar to their non-supersymmetric counter¬ 
parts [15]. Moreover, they are related to the super-Virasoro singular vectors [1, 16]. 

"^By contrast, the AGT-type double Jack polynomials [2, 30] are composed of two Jacks with different variables 
(albeit corresponding to a less radical plethystic transformation), but the same coupling constant. 
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for which these are eigenfunctions^. Somewhat unexpectedly, this Hamiltonian is built in part 
from the generators of the nonnegative modes of an s [2 algebra. 

The article is organized as follows. In Section 2, we briefly review the Calogero-Sutherland 
model and their eigenfunctions, emphasizing a Fock space representation to be used throughout 
the article. The supersymmetric CS model is introduced in Section 3, together with the Jack 
superpolynomials. For a sufficiently high fermionic degree, the supersymmetric Hamiltonian 
eigenvalues are shown to be decomposable into two independent parts. This points toward 
the splitting of the Hamiltonian into two independent CS Hamiltonians and the corresponding 
factorization of its eigenfunctions. The resulting double Jack polynomials are defined formally 
in Section 4 and exemplified for simple cases, while their corresponding Hamiltonian is derived 
in Section 5. 


2 The Calogero-Sutherland model and Jack polynomials 


The CS model describes a system of N identical particles of mass m lying on a circle of cir¬ 
cumference L and interacting pairwise through the inverse of chord distance squared. Setting 
m = h = 1 and L = 27r, the Hamiltonian reads [42, 43, 44]: 


= - 


1 


N 

2 ^ 
i=l 


m - 1 ) 




^ 4sin^ — Xj)' 


where /? is a dimensionless real coupling constant and [xj,pk] = iSjk-^ To the ground state 
correspond the following wavefunction and eigenvalue: 

V’o(a;) = n I sin \{xj - Xk)f with Eq = ^ -31. 

j<k 


It is convenient to define Zj = and to factor out the contribution of the ground state by 
redefining a gauged Hamiltonian as — Eo)'ipo//3 and to set /3 = l/a: 



This is our starting point. 

The symmetric and triangular eigenfunctions of (2.1) are known as the Jack polynomials 
J^\z) [25]^, where the index A stands for a partition A = (Ai, A 2 ,..., Aw), with the Aj’s being 
non-negative integers such that A* > Aj+i. Their eigenvalues are 

= 2an(X') — 2n(A) -|- (N — 1 -|- a)|A|, 

where [27] 


n(A) = ^(i-l)A. = X;(2‘). 

i i ^ ^ 

®We note that the latter aspect would have been very difficult to study for the Macdonald case given the 
complexity of the supersymmetric form of the corresponding Ruisjenaars-Schneider model [10]. 

®See [24] for an extensive and very clear presentation of the CS model. 

^For a physicist introduction to the Jack polynomials, we refer to [21, 24], A more mathematical presentation 
can be found in [27]. 
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Here A' is the conjugate of A obtained from A by replacing rows by columns in its diagrammatic 

representation, and |A| = ^Aj is the degree of A. We will be interested in the behavior of 

i 

the wavefunction when N is large. It is thus preferable to remove the dependency in N in the 
eigenvalue. For this, we note that J^\z) is homogeneous in the Zi’s, so that it is an eigenfunction 
of the momentum operator V: 

i 

Our task is achieved by redefining the Hamiltonian as 

- [N -l + a)'^Zid^^. 


Jack polynomials J^\z) are thus eigenfunctions of with eigenvalues 

4“^ = 2an(A') - 2n(A). (2.3) 

In the large N limit, it is convenient to rewrite the Hamiltonian in terms of power sums 

Pk = Zi+Z 2 -\ -. Since 7^^"^ is a differential operator of order two, it is sufficient to determine 

its action on the product PmPn- A direct computation gives [4, 29] 

77O) = (a - 1) ^ (n^ - n)pndp„ + ^ [(m + n)pmPndp^+„ + amnpm+ndp„dpj. 

n>l n,m>l 

This naturally leads to the Fock space representation 

Q,^(a) = (q, _ 1) _ i)alae + ^ [o^a4fc+£ + > 

e>i k,i>i 


where 


[afc,aj] = ka5k/ and [ak,ai] = [4,4] = 0. 

The correspondence with symmetric functions, together with |0) <—> 1, is 

4 ^—^ Pk Qk <—> kadp^^. (2.4) 

This correspondence preserves the commutation relations. In this representation, the eigenfunc¬ 
tions take the form of a combination of states 

4“^ (4>4>4>---)io)- 

(To be clear, this mode expression is a formal representation of the Jack polynomial expanded 
in the power-sum basis through the correspondence (2.4).) For instance, up to a multiplicative 
constant® 

44i)|0) [(4)^ + (3® “ 1)4(4)^ + 2a(a — 1)44 “ “(4)^ “ 2a^4] |0)- 

As a side remark, we point out that it is through the correspondence (2.4) that the connection 
between Virasoro singular vectors and Jack polynomials is established [4, 5, 28, 34, 36]. 

®In the monic normalization = m\ + lower terms, where m\ is the monomial symmetric function, this 
coefficient is 1/[2(1 -|- a)^]. 
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3 Supersymmetric version 


In order to supersymmetrize the CS model, we need to introduce anticommuting variables 
9i,... ,9 m and extend the CS Hamiltonian H in the following way: 

'H^uly = {Q, + terms depending upon 9i, 

for two fermionic charges Q and of the form Q = Y2 and dQ^A\{x,p), where 

i i 

Ai and A\ are fixed by the requirement of reproducing the 7^^“^ term on the rhs of the above 
equation. This construction leads to 


— 9 

' ^susy ' ^ ^ 


E 




{Zi - ZjY 


{9i - 9j){d8^ - dej). 


This operator is part of the tower of conserved quantities Tin, 1 < n < N {V = T-Li and T^susy = 
7 ^ 2 ) that reduce to the usual (gauged) CS conservation laws in the absence of anticommuting 
variables. But given that there are 2N degrees of freedom in the supersymmetric version, there 
are N extra conserved charges that vanish when all 9i = Q [14]. The first nontrivial representative 
of this second tower is 


t(“) 

-^susy 


N 


= a 


1<2<J<A^ 


2 = 1 


ZjOi 
Zi - Zn 


{de. 


de,). 


As a side remark, we mention that both expressions can be represented in the Fock space of 
a free boson, described by the modes Ufe, a\ (with fc > 1, i.e., without the zero mode) and a free 
fermion, whose modes are denoted bk, 


«^iuiy = («-!) ~ + Y 

r>i fc,r>i 

+ a{a - 1) ^ - ^)blbe + aY ‘^^[alb\bk+(. + bYbiUk] 

£>1 k,i>l 


and 


= (« - 1) X] + X] + akblbi+k] ■ 

l>0 £>0,k>l 

The fermionic modes are governed by the anticommutation relations: 

{bk,bl}=6k/ {bk,b£} = {bl,bj} = 0, 

and their correspondence with symmetric functions is 

bl <—>Pk = 0izY^2zY - and bk <—> dp^^. (3.1) 

Now, assuming a natural triangularity condition, the common eigenfunctions of T^iuiy and zYy 
are the Jack polynomials in superspace, or Jack superpolynomials, denoted by jY(z, 9) [14]. 
They are homogeneous in z and in 9 and invariant under the exchange of pairs {zi,9i) i—> 

®In a supersymmetric context, the modes of the partner free fermion should pertain to the Neveu-Schwarz 
sector, hence be half-integers. This can be achieved by redefining (bt,b],) as (bfe+ 1 / 2 , relation (3.1) 

below. However, this precision is not required in the present context. 
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{zj,9j). Their labelling index A is a superpartition. Before displaying the eigenvalues, some 
notation related to superpartitions is required. 

A superpartition A is a pair of partitions 


A = (A“; A^) such that 


A^ is an ordinary partition, 

A“ is a partition with no repeated parts. 


Note that the last part of A“ is allowed to be zero. We denote by A* the partition obtained by 
reordering in non-increasing order the entries of A“ and A^ concatenated. The diagrammatic 
representation of A is obtained by putting dots at the end of the rows that come from A“ (in 
such a way that dots never lie under an empty cell). Here is an example: 


A = (4,2,0; 3, 2,1,1) ^ (4, 3, 2, 2,1,1,0) ^ 


A superpartition is equally well described by the pair A* and A®, where the latter is the partition 
obtained by replacing dots by boxes, e.g., in the example above. 


A* = 



Finally, the bosonic degree of a superpartition is the number of boxes of A* and the fermionic 
degree, generally denoted by m, is the number of dots in the diagram of A, that is, the number 
of parts of A“. 

We are now in position to give the eigenvalues of 'Hinly and Isuly corresponding to the 
eigenfunction ■ These are respectively 


= 2an{h*') — 2n(A*) and = q;|A“| — |A'“ 


In the supersymmetric case, we are not only interested in the large N limit but also in the 
large m limit (actually, in the large m and N—m limits). We thus want to extract from the above 
two eigenvalues, their dependence on m which is somewhat hidden. For this, we first notice that 
when m is large (relative to the size of A, an estimation that is made precise in (3.4)), there 
are circles in every possible positions in the diagram of A.^*^ As such, the circles can be ignored 
and we observe that A* differs slightly from its core = (m — 1, m — 2,..., 1, 0). In the 

diagrammatic representation of A*, the deviations to the core are located at the top right and at 
the bottom left of the diagram. We thus see that the superpartition can be disentangled into its 
fermionic core plus two small partitions A and /r such that A = (A -|- /r) [9].^^ For instance, 

for m = 8, we have 


A = 




^°Here is a technical precision that could safely be skipped. There are unimportant exceptions to the statement 
that when m is sufficiently large (meaning larger or equal to its lower bound, which is |A| + \fi\ for A and n 
defined below), there are dots in every possible positions. That all allowed slots are filled by dots is true when 
m > £(A) + 1 + j-Li, i{X) being the length of the partition A. Since |A| + |/r| + 1 > + 1 + /Ui ths statement is 

always true for instance when m>|A| + |/i| + l. 

^'^For the + operation, the parts add up. For example, we have (3,1) + (4, 2, 2) = (7, 3, 2). 
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It is clear that A is fully characterized by m and the pair (A, /r) (whose total degree is much 
less than that of A). The main advantage of this diagrammatic decomposition is that it implies 
readily that when m is large the conjugate of A is A' = A'). 

Let us reformulate the eigenvalues in terms of the data A, ^ and m. For the ZsuJy eigenvalue, 
the computation is easy and yields 


(a) m large (a) 


"A 


= "l-^l “ + (« - l)m(m - l)/2. 


We can easily remove the dependency in m in the eigenvalue by redefining xiusy as follows: 


4uiy ^ ^iuiy = ^iuiy ” ” 1)M(M - l)/2, where M = Y^ 0.9,, 

i 

This subtraction is well defined since AJ is also a conserved quantity. The modified eigenvalue 
reads then 


At = “l-'l - M- (3.2) 

The eigenvalues of L^luiy can also be reformulated in terms of A, /r and m, again keeping in 
mind that this is valid only for sufficiently large m. Observe that^^ 

A* = (a + U /r => A* G {Aj + rn — j\l<j< mj U {/i^ | 1 < A: < 


and similarly for A*^ = + S^"^^)L)X', where £(/r) is the length of the partition (the number of 

non-zero parts). Fortunately, the calculation of n{A*) (and n{A*')) is independent of the precise 
relationship between the indices i and j, k in the above notation if we use the second expression 
of n(A) given in (2.2). Let us first consider 


n{A*') 


E 

i 

m 

E 

2=1 


E 


Aj + m — i 


+ 


2=1 

— i 


+ Aj(m 






ii\) 

n(A') +n(/i') -b ^ Xi[{m — 1) — (i — 1)] + m{m — l){m — 2)/6 
2=1 

n(A') -b re(/i') -b (m — 1)|A| — n(A) + m{m — l)(m — 2)/6, 


where in the last step, we use the first expression in (2.2). For the computation of n(A*), we 
simply replace A and fj. by fi' and A' respectively in the previous expression to get: 

n(A*) = n(/i) + n(A) + (m — l)|/r| — n(/i') + m{m — l)(m — 2)/6. 


Combining these two expressions yields 


(a) m large 


= (a + -p 2(m — 1)6^°"^ + (a — l)m(m — l)(m — 2)/3, 


where is defined in (2.3). We can thus remove the dependency in m in the eigenvalue by 
redefining L^iuiy as 

^iuiy ^ ^iuiy = ^iuiy ” 2(M - 1)1^% - (« - l)M{M - 1){M - 2)/3. 

^^For the U operation, the rows of the second partition are inserted into the first one; for instance (3,1) U 
(4,2,2) = (4,3, 2, 2,1). 
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The ?iiusy eigenvalue is then simply 

4'2 = (“ + l)4“''“*"'+e<r+‘'. (3.3) 

The two expressions (3.2) and (3.3) are the results we were looking for, the conclusion of our 
heuristic argument based on simple transformations of the two defining Jack-superpolynomials’ 
eigen-operators. On the one hand, that the m-dependence of the eigenvalues can be removed is 
an indication of the stability property of the eigenfunctions. On the other hand, the decoupling of 
these two eigenvalues into two independent sectors A and /i suggests the factorization property. 
More precisely, the decomposition of as the sum of and hints that, in 

the large m limit, the eigenfunction should somehow factorize into a product of the form 

(with (3.2) being compatible with the respective modifications of the 

coupling constants). 

These two conclusions are indeed confirmed: the eigenfunctions both stabilize and factorize 
(after a certain transformation that will be explained in equation (3.5))^^ for [9]: 


m>\X\ + \fj.\ 


(3.4) 


Let us consider a simple example. For {X, fi) = ((!),(!)) the m = 1,2, 3,4 eigenfunctions 
read respectively 

[^o(«i)^ + Oib\a\ - blal - abl] |0), 

[b\bl{a\f + {a- l)blbla\ - ablb\ - ablbl] | 0 ), 

[blb\bl{a\)‘^ + {a- l)blb\bla\ - - ablb\bl] |0), 

-b (a — l)b\b\b\bQa\ — ab4b\b\bQ — 06562^^61!,] |0). 


Clearly, the m = 1 (< |A| -|- |/i| = 2) wavefunction does not belong to the stable sector. For 
m> 2, the coefficients and the content of each term are always the same. This is the stability 
property. 

Although they stabilize, the eigenfunctions still depend on m. However, consider the map 
for (A“; A^) ^^ (A,/x): ^ sxiy)Pf,{y, z), (3.5) 


where 

m N—m 

PtJ.e with pn{y, z) = J2yi + 

i=l i=l 

and s\ is the Schur function. Observe that Pn{y, z) is simply pn in the variables yi,y2, ■ ■ ■, ym, zi, 
Z 2 ,... ,Z]sf-m- This maps the above eigenfunctions corresponding to the values m = 2,3,4 to 
(inserting the proper normalization) 

^ (1 + a) + {a- l)si{y)pi{y,z) - asi^i{y) - as2{y)]. 

The stability has now been lifted to the full structure of the eigenfunction. 

When the Jack superpolynomial is expressed in terms of the variables {x, 9) rather than in modes, the 
transformation is simply 

Amix)~^de^ ■ ■ ■ deiJj^\x,6), where A^(x) = ](([ (xi-Xj). 

l<2<j<m 
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But in addition, the map (3.5) captures the factorization property suggested by the form of 
the eigenvalues. Using the Fieri rule for Schur functions [6, 27] to express the sum of the last 
two terms in a product form, 

si,i{y) + S2{y) = si{y)si{y), 

we see that can also be written in a product form 

= JY^r^{Piiy^z) + asi{y)){pi{y,z)-siiy)). (3.6) 

This is a simple illustration of the announced factorization. 


4 Double Jack polynomials 


In general, the action of the map (3.5) at the level of Jack polynomials is [9] 

where 




(q/(q+1)) 


y +-z 

a + 1 


j(“+i)( 


z ■ 


(4.1) 


Here we use the plethystic notation (see, e.g., [6, 26]). In our case, it simply means that if 
J_[“/(“"'~i))(pi,P 2 ,P 3 , •. •) is the expression of in terms of power-sums, then 


J 


(“/(«+!)) 


y + 


ol \ 


-Z 


= (piiy) + —^pi{z),p2{y) + —^P2{z),Pz{y) + —^P2.{z),.. ■ 

\ a + 1 a + \ a + \ 


that is, z/{I + a)] is obtained from the expansion of in terms of 

power-sums by replacing by pn{y) + -^Pn{z). 

Let us recover (3.6) from the general expression (4.1). This is a particularly simple case given 
that = S(i) = 771 ( 1 ) = Pi- With pi{y, z) = pi{y) +pi{z), J(")^(^)(y, z) becomes 

= (pi{y) + 


which is indeed of the form (4.1). Here is a slightly more complicated example: 

r(«) ' ' ^ 




{ 2 ),(o)v«>~/ (i + Q,)(i + 2a) 

With S( 2 ) = {pi +p 2 )/‘i, simple algebra yields 


[pi{y, zf + ap2{y, z) + 2 as(i)(y)pi(y, z) + 2a^S(^2){y)] ■ 


4"uoi(yw) = 


(1 J- a) 


a 


kW + ) = L) 


(«/(«+!)) 


[^] 


'{2),{0)V«>~/ (l + 2a) 

with X = Y+Z/ (a-Fl) and where in the last step, we used the expression = {Pi+Cip 2 ) /(U-a). 
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A more formal characterization of z), which we call the double Jack polynomials, is as 


follows [9]. They are the unique bi-symmetric functions such that 


z) = s\{y)s^{z) + smaller terms 


and 


= 0 if (A,^) / (z^,k). 


The triangularity condition that specifies the “smaller terms” refers to the double version of the 
dominance ordering: 


{X, fj.) > {n, k) ^ 

e 

- ni) > 0, 


|A| + \y\ — \i^\ + |k|, 


— \1^\ 


2 = 1 


i=i 


VT 


while the orthogonality condition refers to the scalar product 
{sx{y)PiM{y, z), Su{y)pi^{y, z)) = 


(4.2) 


with z^= W n^(i) being the multiplicity of the part i in p. Observe that this scalar 

i>l 

product has the form 




Jack’ 


where (•)')gj,hui. is the scalar product with respect to which the Schur functions s\{y) are or¬ 
thonormal while (•, is the scalar product with respect to which the Jack polynomials 

J^\y, z) are orthogonal {J^\y, z) being the usual Jack polynomials in the variables yi,y 2 , ■ ■ ■, 

Um 1 ^2-1 —m) ■ 


5 The double CS model and an emerging SI 2 

Let us now unravel the integrable model whose eigenfunctions are the double Jack polynomials. 
The factorized expression (4.1) of these polynomials and the splitting of the eigenvalue displayed 
in (3.3) readily indicate that the underlying Hamiltonian %d a. sum of two CS Hamiltonians, 
albeit with modified coupling constants and involving unusual variables: 


JId = (a + l)J2i + 


where 


Hi = hW(“+i)) with 


Pn ^ Pn[X], 


9p„ ^ d. 


where X = Y + {a + 1) ^Z, 


Pn [^] ’ 


and 


H 2 = 


with 


Pn ^ Pniz), 

^dpn ^ dpn{z)- 


Note that Pn[X] and Pn{z) are considered to be independent. Being the sum of two independent 
integrable Hamiltonians, H_d trivially characterizes a new integrable model. 
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However, the above splitting of is not very interesting since it is hard to give a physical 
meaning to the power-sums Pn[X], Pn{z) and their derivatives. The structure of the scalar 
product (4.2) points toward a more interesting choice of variables, namely Pn{y) and pn{y,z), 
whose adjoints are ndp^(^y^ and nadp^i^y^z) respectively. With 

X = y + = -^Y + —^{Y + Z), 

a + 1 a + l a-br 

the change of variables is thus 

cx 1 

Pn[X] = - —Pn{y) H- ^Pn{y, z), Pn{z) = Pn{y, z) - Pn{y), 

a+l a+l 

which gives (using the chain rule in two variables) 

a 


^Pn[X] ^Pn(y) 


a 


Pn+) 


d, 


1 


These expressions are readily checked by verifying that they satisfy the commutation relations: 


-d, 


Pn(y)- 


(5.1) 


(5.2) 


J — ^n,mj [^Pn(2)) Pm(^)] — ^n,rm 

[dp+X],Pm{z)] = 0, [dp^^^),Pm[X]] = 0. 

For these manipulations, we stress that Pn{y) and pn{y,z) are considered to be independent, 
meaning: 


[dp„(y),Pm{y,z)] = [dp^(y,z),Pm{y)] = 0 . 

We then substitute (5.1) and (5.2) into (a+l)?^i+?^ 2 - The result, obtained after straightforward 
manipulations, is best rewritten in terms of two independent sets of bosonic modes defined as 

^l=Pn{y) and An = ndp^(^y) (^ [Hfc, hJ] = Mfc,r), 

together with 

ah=Pn{y-,z) and Un = nadp^^y^^) (^ [ak,al] = ka6k/) ■ 

The resulting form of T-Ld is 

a+Ln = ^ [a\^\ak+i + + (a - 1) ^(^ - l)ajar - a ^(^ - 1) + Hjar] 

fc,r>i r>i £>i 

+ oi \2a}^A^f^Ak+i + a}^j_f^AkA(^ + a + 2j4j,_|_^Tfcaf] 

k,i>i fc,r>i 

+ a(a-l)^ [4 A^iAk-\-(, ^k+e^kAt\. (5.3) 

fc,r>i 

It turns out that "Hd can be reexpressed as 

Hd = +“> + (a - l)+‘> + [Qi.w!,"] - \[QuQ2], 


ei = E)[“+(- 


and 


r>i 


Q2 = Y,{e-l)[AlAe-^al 


e>i 


where 
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Note that 


IS 


in the variables yi, 


2 / 2 , 


, 2/m, ■2-1, -22, • • • , Zpf—mi SO that 


^ [alalttk+i + al^^Ukae] + {a - 
k,e>i 


1 )^(£- l)a\at. 
£>l 


Similarly, is in the variables 2/i, 2/2, • • ■, 2 /m but evaluated at a = 1: 
a(a - = a{a - 1) ^ [-^l^l^k+e + Al_^_gAkAi]. 


kl>l 


Next, it is simple to check that yields the second line in (5.3). Therefore, parts 

of the constituents of T-Ld have a direct interpretation in terms of variables. However, this is 
not the case for Qi and Q 2 - Note that the action of Qi amounts to exchanging the a and A 
modes (which thereby appears to be a remnant of the action of a supersymmetric charge). 
Nevertheless, it turns out that Qi and Q 2 have a nice Lie algebraic interpretation. More precisely, 
both are combinations of the generators of an underlying affine sb algebra (whose existence 
is not surprising in the presence of two independent infinite sets of bosonic modes). It is 
straightforward to verify that the operators 

e>i £>i £>i 

do satisfy the sb commutation relations 

[eW,/W] = [/.W,eW] = = -2/("+^). 

We thus get that 

Qi = and Q 2 = 

and, as such, T-Ln Is built from a special intertwining of i-L^} and i-6y'^ with the genera¬ 
tors and of the nonnegative part of 5b- 

This intertwining pattern is expected to hold for all the conserved quantities of the double 
CS model. Consider for instance the two conserved quantities of degree 1 


A^^Ai - a\ai 


Id — (cK ~ 1) Al^A i + [ajA£ + Aja^J, 

i>l £>1 


Id = Y 


£>1 


A^ A£ -I — a£ 


a 


whose eigenvalues are respectively (3.2) and |A| -|- \fi\. As for T-Id, the conserved quantity Id can 
be written in terms of the usual conserved quantities of the two CS models specified by 
and and the (nonnegative-mode) generators of sb: 

Id = {d — f)'Py + [Qli ly ], 

where Vy is the momentum operator V in the variables yi,... ,ym- Similarly, we have 


Id — ly + Iy,z- 
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A The many combinations of “super”, “symmetric” 
and “polynomials” 

The supersymmetric CS model involves, in addition to the usual N variables Xi, their anti¬ 
commuting partners As a result, its eigenfunctions depend upon both set of variables. To 
describe a function with such a dependence on both types of variables, we use the terminology 
“superfunction”. The qualifier “symmetric”, when applied to superpolynomials refers to the 
invariance upon the interchange of the pairs {xi,9i) and {xj,9j), an entanglement which is in 
the spirit of supersymmetry, as opposed to invariance under Xi -H- Xj and 9^ -H- 9i separately. In 
this sense, it is thus natural to call the model’s eigenfunctions the Jack superpolynomials. 

Unfortunately, this terminology is a bit ambiguous with regard to existing similar nomencla¬ 
tures. For instance, the symmetric superpolynomials just defined are completely different from 
the “supersymmetric polynomials” considered, e.g., in [41] (which incidentally have nothing to 
do with the physicists’ concept of supersymmetry and do not involve anticommuting variables). 
Indeed, the latter are defined by two conditions: (1) They are doubly symmetric polynomials 
in two distinct sets of ordinary (commuting) variables xi,, Xm and yi,..., i.e., invariant 

under independent permutations of the Xj’s and the y^’s. (2) They satisfy the following cance¬ 
lation condition: by substituting xi = t and yi = t, the polynomials become independent of t 
(this is the meaning of “supersymmetric” in this context). An example of a generating function 
for such polynomials is 

m n 

]J(1 - qxi) ]J(1 - qyj)~^ = '^P{r){x,y)q''■ 

i=l j=l r>0 

The terminology “super-Jack polynomials” is used by the authors of [37, 38, 39] referring 
to objects also different from our Jack superpolynomials. Their construction originated from 
the realization that the CS model has an underlying Ajv-type root structure. Viewed from 
this perspective, the model can be generalized to arbitrary root systems by still preserving 
integrability [32, 33] . The models considered in [38, 39] are natural extensions of this construction 
but based on Lie superalgebras (see also [12]). Here again, there are no anticommuting variables. 
By contrast, our Jack superpolynomials have no known Lie-algebraic interpretation. 

Still another supersymmetric generalization of the CS model is presented in [22, 23], albeit not 
formulated in terms of a supersymmetric quantum mechanical problem. Here the Hamiltonian 
is induced from the radial part of the Laplace-Beltrami operator on symmetric superspaces. 
However, it does not contain anticommuting variables. The resulting models are a one-parameter 
generalizations of those of [38, 39], but with the integrability issue unsettled. 

To complete this list, we note that the term “superpolynomial” is used in the theory of knots, 
with a completely different meaning (it refers to the Poincare polynomial for a quite general - 
triply-graded - homology) [18]. 

To avoid terminological confusion, we generally employ the nomenclature “Jack polynomials 
in superspace”, using “Jack superpolynomials” as an abbreviation. But where is the superspace? 
Indeed, the usual variables Xi do not refer to ordinary space. However, the Jack polynomials 
are wave functions of a V-body quantum mechanical problem and as such they are defined in 
the configuration space of the various particles described by the positions x,, i = 1,..., V. It is 
this configuration space that is lifted to superspace. 
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